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1.0  INTRODUCTION 


This  report,  prepared  by  the  Grumman  Corporate  Research  Center  under 
Contract  No.  DAAE07-86-C-R020,  deals  with  the  dynamic  response  of  soil 
to  loads  generated  by  the  tracks  of  military  vehicles.  The  simulation 
of  the  dynamic  behavior  of  vehicles  traveling  off-road  is  of  great 
interest  to  the  Army,  since  vibrations  resulting  from  the  unevenness  of 
the  terrain  adversely  affect  the  ability  of  the  crew  to  operate  the 
vehicle  and  its  weaponry.  Increased  battlefield  mobility,  coupled  with 
the  demand  that  tanks  be  able  to  fire  on  the  move,  requires  gun  control 
systems  which  are  stabilized  irrespective  of  the  motions  of  the 
vehicle.  A  capability  to  simulate  the  dynamic  behavior  of  vehicles  and 
its  components  is  essential  for  the  design  of  these  vehicles  and 
optimization  of  vehicle  components  to  reduce  harmful  vibrations. 


2.0  OBJECTIVES 

The  objectives  of  this  work  are  to: 

•  Evaluate  the  experimental  information  on  dynamic  soil 
response  to  plate  penetration  with  respect  to  its 
relevance  to  dynamic  track-soil  interaction 

•  Determine  the  applicability  of  state-of-the  art  methods  of 
soil  dynamics  to  the  modeling  of  dynamic  soil  response  to 
track  segment  loads 

•  Develop  a  concept  of  dynamic  track-soil  interaction,  to  be 
used  in  conjunction  with  tracked  vehicle  dynamic  models  for 
the  simulation  of  this  interaction  by  computer  models 

•  Write  a  computer  program  as  a  subroutine  to  the  dynamic  track 
model  under  development  at  TACOM  for  the  calculation  of  the 
dynamic  soil  response  to  track  segment  loads. 


3.0  CONCLUSIONS 

The  dynamic  response  of  soils  to  plate  loads  has  been  examined  and 
analyzed  in  detail  to  establish  the  significant  factors  in  the  complex 
interaction  between  track  segment  and  dynamic  soil  response.  These  were 
found  as  follows: 

•  The  loading  history  of  the  soil 

t  The  time  rate  of  load  application 

•  Tractive  forces  developing  shear  stresses  at  the  track-soil 
i nterf ace 
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•  Interaction  between  adjacent  track  segments. 

A  concept  of  dynamic  soil  response  model  has  been  formulated  which  takes 
into  account  all  significant  factors  enumerated  above.  A  computer  code 
for  a  subroutine  suitable  to  use  with  dynamic  track  models  was 
written.  The  limitations  on  CPU  time  allocated  for  the  computation  of 
dynamic  soil  response  in  the  dynamic  track  model  were  met  by  performing 
•all  time-consuming  computations  in  a  preprocessor  program  and 
establishing  in  that  program  appropriate  approximations  for  the 
interrelationships  affecting  the  soil  response. 


4.0  RECOMMENDATIONS 

4.1  Basic  Research 

It  is  recommended  that  a  long-range  research  project  be  initiated  and 
supported  to  advance  the  understanding  of  the  dynamics  of  partially 
saturated  surface  soils  and  develop  theories  to  describe  soil  response 
to  dynamic  loadings  by  off-road  vehicles  in  terms  of  fundamental  soil 
properties. 

4.2  Dynamic  Plate  Sinkage  Tests 

It  is  recommended  that  laboratory  and  field  testing  apparatuses,  capable 
of  performing  plate  sinkage  tests  in  soils  at  various  penetration 
velocities,  be  designed  and  constructed  and  that  dynamic  plate  sinkage 
tests  be  made  a  standard  procedure  in  all  field  tests  involving  off-road 
vehicles.  A  data  base,  using  the  results  of  these  tests,  should  be 
established  for  the  estimation  of  soil  dynamic  parameters  from  other 
terrain  descriptors. 


5.0  DISCUSSION 

5.1  Review  of  Theories  and  Experimental  Information 

Theoretical  treatment  of  the  behavior  of  soils  under  dynamic  loading 
conditions  requires  the  establishment  of  constitutive  relations  for 
nonlinear  stress-strain  properties  which  depend  both  on  the  loading 
history  and  the  time  rates  of  loading.  Of  the  few  constitutive 
relations  proposed  to  simulate  nonlinear  stress-strain  behavior 
experienced  in  trlaxial  tests,  none  meets  the  test  of  generality, 
although  some,  applied  to  specific  type  of  soil  and  loading  conditions, 
represent  the  behavior  of  that  type  of  soil  fairly  well.  Even  if  a 
general  constitutive  relation  were  available,  its  implementation  into  a 
nonlinear,  three-dimensional  finite  element  program  would  require  a 
major  research  effort,  clearly  outside  of  the  scope  of  the  present  work. 

Theories  of  dynamic  soil  behavior  with  a  limited  scope  have  been 
developed  in  connection  with  the  following  engineering  fields: 
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•  Effects  of  earthquakes  on  structures 

•  Machine  foundations 

•  Pile  driving  dynamics. 

While  the  theories  and  experimental  information  developed  in  connection 
with  these  problem  areas  are  relevant  to  the  problem  of  track-terrain 
interaction,  they  are  not  directly  applicable  to  it,  primarily  because 
they  refer  to  the  dynamic  response  of  soil  at  some  depth  as  opposed  to 
the  surface  loading  of  soil  by  vehicles.  Also,  the  experiments  on  soil 
behavior  performed  in  connection  with  these  problem  areas  are  generally 
of  the  cyclic  rather  than  transient  impulse-type  loading  and  have, 
therefore,  no  direct  application  to  the  problem  of  soil  response  to 
track  action. 

Experiments  performed  at  the  Waterways  Experiment  Station  (WES)  in 
connection  with  off-road  mobility  research  and  at  various  other 
institutions  under  the  sponsorship  of  the  Air  Force  in  connection  with 
aircraft  landing  on  unprepared  airfields  have  direct  bearing  on  the 
problem  of  dynamic  track-soil  interaction.  Table  5-1  summarizes  the 
test  programs  which  were  found  relevant  and  useful  for  the  problem  of 
dynamic  track-soil  interaction.  These  were  analyzed  in  detail;  results 
of  the  analyses  are  presented  in  the  relevant  portions  of  the  subsequent 
discussions. 

5.2.  Principal  Considerations  in  the  Conceptual  Development  of  Dynamic 
Track-Soil  Interaction  Model 

In  the  conceptual  development  of  the  dynamic  track-soil  interaction 
model,  it  was  necessary  to  reconcile  the  requirements  of  realistic 
simulation  with  the  limitation  imposed  by  the  CPU  time  allocated  in  the 
dynamic  track  model  to  soil  response  computations.  Finite  element 
modeling  (FEM)  of  the  dynamic  soil  response  would  be  desirable  from  the 
viewpoint  of  realistic  simulation1,  yet  it  is  prohibitive  because  of  the 
considerable  CPU  time  used  by  FEMs.  In  addition,  FEMs  generally  solve 
for  the  inverse  of  the  dynamic  problem,  namely,  they  calculate  the 
displacements  of  soil  for  applied  forces  as  input. 

To  resolve  the  problem  of  realistic  simulation  within  the  allocated  CPU 
time,  the  following  approach  was  adopted.  The  basic  element  in  the  soil 
response  model  is  the  pressure-sinkage  relationship  obtained  when 
loading  the  soil  by  a  rectangular  plate  of  the  size  of  a  track 
segment.  This  relationship  may  be  represented  by  a  mathematical 
expression  which  lends  itself  readily  to  the  solution  of  the  problem  as 
formulated  in  the  dynamic  track-model:  calculation  of  soil  response 
forces  for  given  track  segment  displacements.  Plate  sinkage  tests  at 
low  penetration  rates  have  been  regularly  conducted,  both  in  the 
laboratory  and  in  the  field,  for  off-road  mobility  research  and 
analyses.  Therefore,  pressure-sinkage  curves  may  easily  be  determined, 
and  the  parameters  in  the  mathematical  expressions  may  be  easily 
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Table  5-1.  Dynamic  Plate  Penetration  Tests 


Institutions: 

WES 

I  IT 

U.  Of  DAYTON 

Experiments  in  sand 

Properties: 

Yuma  Sand 

Newport  News  Riverwash 

sand 

sand 

Median  diameter  (mm) 

0.12 

N.A. 

0.8 

Moisture  content  (%) 

Dry 

8.1 

4.75-5.9 

Dry  unit  weight  (pcf) 

92-107 

99.8** 

113-125 

Degree  of  saturation*^) 

0 

32 

31-37 

Cohesion*  (psf) 

0 

10 

17 

Friction  angle 

39°-46° 

?N.A 

N.A 

Experiments  in  clay 

Alluvial  Mississippi  (Buckshot)  clay 


Moisture  content  (%) 

29.2-39.8 

30.4 

29.4 

Wet  density  (pcf) 

111-118 

105.5** 

107.8** 

Dry  density  (pcf) 

79.5-91.4 

__  _** 
80.7 

__  .** 

83.4 

Degree  of  saturation  (%) 

95 

90 

90 

Cone  index  (psi) 

30-116 

N.A. 

120** 

* 


estimated 

average 
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evaluated.  The  dynamic  effect  at  higher  penetration  rates,  the  effect 
of  tractive  forces  on  the  pressure-si nkage  relationship,  and  the  effect 
of  adjacent  track  segments  may  be  considered  by  appropriate  factors,  or 
modifications  of  the  parameters  in  the  pressure-sinkage  equation.  An 
important  consideration  in  the  conceptual  development  of  the  dynamic 
soil  response  model  has  been  to  take  into  account  the  effect  of  the 
loading  history  of  soil  on  its  dynamic  response.  Soil  is  essentially  a 
material  with  memory;  it  "remembers"  the  stresses  and  displacements 
imposed  on  it  by  previous  actions  when  responding  to  a  new  excitation. 
Since  the  soil  is  loaded  and  unloaded  again  and  again  during  the  pass  of 
a  track,  it  is  essential  to  consider  the  changes  in  soil  response  as  the 
various  track  segments  pass  over  it. 

5.3  Pressure-Sinkaqe  Relationship  for  Individual  Track  Segments 

Pressure-sinkage  relationships  have  been  used  in  off-road  mobility 
research  to  characterize  soil  behavior  since  the  late  fifties.  The 
exponential  relationship 

p  =  k.zn  [1] 

where  p  =  pressure 
z  =  si nkage 
k,n  =  parameters 

and  its  modifications,  advocated  by  Bekker  for  the  determination  of  the 
pressure  distribution  beneath  wheels,  characterize  the  soil  behavior  by 
the  two  parameters  k  and  n.  Since  the  dimension  of  the  parameter  "k" 
depends  on  the  exponent  "n,"  it  is  difficult  to  associate  a  "k"  value 
with  a  physical  property  of  the  soil,  or  compare  "k"  values  obtained  in 
various  soils.  • 

A  hyperbolic  type  of  equation  proposed  by  Kondner  and  Krizek^  reads  as 
follows: 


p  =  q  - _  12] 

z  (B  +  R  z//A 

where  q^  =  tri axial  (or  reference)  strength  of  soil 
A  =  area  of  plate 
B,R  =  parameters 

In  this  equation,  the  dimensionless  coefficients  B  and  R  characterize 
the  behavior  of  the  soil  together  with  the  value  of  q^,  the  triaxial 
strength  of  the  soil,  or  a  selected  reference  strength.  The  coefficient 
B  is  proportional  to  the  tangent  of  the  pressure-sinkage  curve  at  z  =  0, 
the  coefficient  of  proportional ity  being  1/q^. .  The  coefficient  R  is 
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defined  as  follows: 


B 

z//A 


[3] 


R  = 


Experience  with  plate  penetration  tests  has  shown  that  at  penetrations 
in  the  range  of  z//A  >  0.5,  the  soil  is  in  the  plastic  state.  The 
second  term  in  Equation  [3]  then  becomes  very  small,  and  the  value  of  R 
may  be  associated  with  the  approximate  value  of  the  ratio  of  the 
triaxial  strength  to  the  bearing  strength. 

To  determine  the  applicability  of  Equation  [2]  to  the  pressure-sinkage 
curves  obtained  in  dynamic  plate  penetration  tests,  the  experiments 
listed  in  Table  5-1  were  analyzed  in  detail.  The  main  findings  of  these 
analyses  are  summarized  hereafter. 

The  graphically  reported  results  of  74  pressure-controlled  plate  sinkage 
tests,  done  at  I IT  Research  Institute3,  were  digitized  and  analyzed  by 
various  computer  programs  to  determine  the  degree  of  correlation  with 
various  parameters  Influencing  the  dynamic  soil  response.  For  example, 
correlations  between  the  developed  pressure  and  rate  of  penetration  were 
found  to  be  very  poor,  indicating  that  linear  viscosity  does  not 
directly  apply  to  the  problem.  In  another  analysis  it  was  found  that 
while  correlation  between  the  sum  of  pressure  times  the  duration  of 
loading  was  good  within  a  certain  group  of  tests,  others  negated  this 
correlation. 

This  test  series  was  also  analyzed  to  determine  how  well  the  proposed 
pressure-sinkage  equation  fits  the  experimental  data.  Equation  [3]  may 
be  rearranged  and  a  new  variable 


p./A 

may  be  introduced,  resulting  in  the  linear  equation  as  follows: 

B  +  ^  =  Y  [5] 

✓A 

Figures  5-1  and  5-2  show  typical  data  points  from  this  test  series 
plotted  against  the  new  variable  Y.  Clearly,  the  data  points  fit  a 
straight  line  fairly  well  for  the  tests  performed  in  both  sand  and  clay, 
except  for  a  few  initial  points  in  Figure  5-1.  These  are  believed  to  be 
the  result  of  seating  inaccuracies  rather  than  signifying  a  different 
behavior  at  very  small  sinkages. 

Equation  [2]  also  fits  reasonably  well  the  results  of  penetration- 
controlled  tests.  Figure  5-3  shows  data  points  from  the  tests  done  at 
the  University  of  Dayton  in  clay4  (see  Table  5-1).  The  curve  fits  by 
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Figure  5-3.  Penetration-controlled  Plate  Sinkage  Tests  on  Clay 
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Figure  5-4.  Results  and  Evaluation  of  Experiments  Performed  at  WES  on  Yung  Sand  (Plate  Size:  1x4  in.) 


16 


"  SCALE  x  lO-*2  (1/PSI) 


Equation  [2]  are  shown  by  dashed  lines.  While  the  results  of  most  of 
the  plate  sinkage  tests  may  be  represented  by  Equation  [2]  fairly  well, 
some  tests  performed  in  sand  show  the  phenomenon  of  brittle  failure, 
caused  by  the  strain  softening  stress-strain  behavior  of  dense  sands. 
Figure  5-4  shows  the  results  of  penetration-controlled  experiments  from 
the  test  series  done  at  WES5.  Dotted  lines  show  an  approximation  of  the 
upper  portion  of  the  observed  pressure-sinkage  relations  by  Equation 
[2].  While  the  approximation  in  the  range  of  brittle  failure,  where 
pressures  drop  with  increasing  sinkage,  is  poor,  it  would  be  necessarily 
so  with  any  other  smooth  monotonic  mathematical  relationship  which  could 
be  considered  for  representation  of  the  soil  response  in  the  dynamic 
track  model.  Since  surface  sands  encountered  in  off-road  travel  are 
highly  unlikely  to  be  in  a  dense  state  comparable  to  the  80%  relative 
density  to  which  the  sand  in  the  WES  tests  shown  in  Figure  5-4  was 
compacted,  it  was  not  deemed  necessary  to  use  a  different  type  of 
pressure-sinkage  relationship  of  great  complexity  for  the  better 
representation  of  this  special  case. 

Shear  stresses  developed  by  tractive  or  braking  forces  at  the  track-soil 
interface  interact  with  the  soil  response  and  modify  the  pressure 
sinkage  relationship  represented  by  Equation  [2],  In  the  plastic  state 
of  soil  this  effect  manifests  itself  in  the  reduction  of  bearing 
stresses.  It  was  found  expedient  to  express  the  effect  of  shear 
stresses  in  terms  of  the  interface  friction  angle,  6,  which  is  the  angle 
enclosed  by  the  resultant  of  the  shear  and  normal  stresses  and  the 
normal  to  the  interface.  The  reduction  factor  F$,  representing  the 
ratio  of  bearing  stresses  at  a  given  6  angle  to  that  at  6  =  0,  may  then 
be  closely  approximated  by  a  parabolic  relationship.  While  this 
reduction  factor  strictly  applies  only  when  the  soil  is  in  the  plastic 
state,  it  has  been  assumed  that  it  approximately  expresses  the  pressure 
reduction  approximately  in  the  elastic-plastic  state  of  soil  as  well. 

The  effect  of  adjacent  track  segments  on  the  pressure-sinkage 
relationship  has  been  analyzed  by  various  methods.  It  has  been  found 
that  even  if  the  interface  friction  angle  is  small,  the  pressure-sinkage 
relationship  is  significantly  affected  only  by  that  track  segment  which 
abuts  a  segment  in  the  direction  of  the  applied  shear  stresses,  and  only 
if  that  track  segment  moves  downward  and  is  in  contact  with  the  ground 
surface  at  that  point.  Under  these  conditions  the  effect  of  an  adjacent 
track  segment  may  be  approximated  by  a  factor  Fa,  which  expresses  the 
ratio  of  bearing  capacity  of  a  plate  twice  the  width  of  a  track  segment 
to  that  of  a  track  segment  standing  alone.  The  dynamic  plate  sinkage 
tests  listed  in  Table  5-1  were  also  analyzed  to  determine  the  effect  of 
penetration  velocity  on  the  soil  response  and  establish  an  approximation 
for  the  representation  of  the  dynamic  effect  which  correlates  well  with 
the  experimental  data.  Of  the  several  concepts  which  were  evaluated, 
the  following  correlated  best  with  the  test  results.  The  "B"  and  "R" 
parameters  in  Equation  [2]  were  evaluated  for  the  plate  sinkage  tests 
performed  at  various  penetration  rates.  It  was  found  that  these 
parameters,  when  plotted  against  the  logarithm  of  the  rate  of 
penetration,  fit  a  straight  line  reasonably  well.  Thus,  the  parameters 
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"B"  and  "R"  were  replaced  by  the  rate  dependent  parameters  "BDYN"  and 
"RDYN,"  and  an  approximate  linear  relationship  was  introduced  for  the 
quick  calculation  of  each  of  these  parameters  in  the  dynamic  track 
model.  These  relationships  are  as  follows: 

BDYN  =  B  +  Fb(Jlog  (z)  [6] 

RDYN  =  R  +  Ffdlog  (z) 


The  coefficients  Fbd  and  Frd  may  be  calculated  by  linear  regression  from 
test  data.  Figures  5-5  and  5-6  show  BDYN  and  RDYN  values,  respectively, 
evaluated  for  the  penetration  controlled  tests  done  at  the  University  of 
Dayton  on  clay  (see  Table  5-1)  and  plotted  against  the  logarithm  of  the 
rate  of  penetration.  The  slope  of  the  best  fitting  straight  line  yields 
the  coefficients  Fbd  and  Frd. 


Track-soil  Interface  for  Dynamic  Track  Model 


5.4.1.  Dynamic  Ground  Response  Plate  (DYGREP).  The  basic  element  of 
the  soil  model  is  the  Dynamic  Ground  Response  Plate  (DYGREP).  For 
computational  efficiency  it  is  expedient  to  use  plates  of  the  same  size 
as  that  of  a  track  segment.  The  position  of  a  plate  is  defined  by  the 
ground  coordinates  of  its  center.  The  DYGREP  plate  can  move  only  in  the 
direction  of  the  z'  coordinate.  (z‘  is  the  downward  coordinate  In  the 
local  plate  coordinate  system  In  accordance  with  ISTVS  nomenclature.) 

The  plate  remains  horizontal  Irrespective  of  changes  in  its  vertical 
position  or  the  elevation  of  adjacent  plates. 


Plate  response  is  activated  by  a  forced  downward  displacement  to  a  new 
position  (z 1 )  of  the  plate.  This  displacement  is  resisted  by  a  uniform 
pressure  over  the  area  of  the  plate.  As  long  as  the  z'  displacement  is 
greater  than  any  previous  downward  displacement  (z'_ ax)  of  the  plate 
(Figure  5-7),  the  magnitude  of  this  pressure  is  controlled  by  the 
pressure-sinkage  relationship  expressed  by  Equation  [2].  The  first  time 
-the  plate  is  forced  downward  from  its  initial  position  (zm,x  =  0), 
Equation  [2]  controls  the  resisting  pressure,  since  z'  >  z?max.  Then, 
before  a  new  z'  value  Is  computed  for  subsequent  time  increments,  the 
value  of  z'max  is  updated  to  z'.  If  there  is  no  displacement  constraint 
on  the  plate,  it  is  assumed  to  rebound  from  its  maximum  displacement 
z'max  as  fon°ws: 

z'o  =  2 'max  -  Fr‘z'max  =  ^  “  Fr>  2'max  [7> 

where  z'0  =  plate  displacement  after  unloading 

z'max  =  max  vertical  displacement  of  the  plate 
Fr  =  coefficient  of  rebound 

If  a  new  displacement  constraint  z'  is  such  that 


z'  <  z' 
0 


<  z 


max 
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Figure  5-5.  Coefficient  "B"  Evaluated  from  Penetration-controlled  Plate  Sinkage  Tests  on  Clay 
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,  Figure  5-6.  Coefficient  "R"  Evaluated  from  Penetration-controlled  Plate  Sinkage  Test*  on  Clay 
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Figure  5-7.  Displacement  —  Pressure  Relationship  for  DYGREP  Plates 
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then  the  resisting  pressure  is  controlled  by  the  relationship  described 
by  Equation  [2],  entered  with  an  adjusted  value  of  z'  as  follows: 

z'adj  ■  <z'  -  z'o>'Fr  I«1 

A  schematic  illustration  of  the  DV6REP  plate  is  shown  in  Figure  5-8. 

The  plate  itself  is  supported  by  nonlinear  springs  and  a  center  post 
connected  to  firm  ground  by  a  ratchet  mechanism.  The  downward  movement 
of  the  plate  is  resisted  by  the  nonlinear  spring  forces,  while  the 
ratchet  mechanism  simulates  the  effect  of  loading  history  on  the  dynamic 
soils  response. 

5.4.2.  Effect  of  Interface  Shear  Stresses  on  Soil  Response.  The  effect 
of  interface  shear  stresses  on  the  soil  response  is  taken  into  account 
by  the  factor  F$,  which  is  approximated  by  the  following  equation: 

Fs  =  FSC(l)  +  FSC(2)  x  Dr  +  FSC(3)  xDf2 


where  FSC(i)  =  coefficients,  i  =  1-3 

Dr  =  interface  friction  angle  development  ratio 

The  interface  friction  angle,  «,  is  the  angle  enclosed  by  the  resultant 
stress  acting  on  the  interface  and  the  vertical.  The  ratio  Dr  is 
defined  as  follows: 


where  6max  =  the  maximum  value  of  the  interface  friction  angle 

The  values  of  the  FSC(i)  coefficients  are  determined  in  the  preprocessor 
program  (Section  5.0). 

5.4.3.  Effect  of  Adjacent  Track  Segments  on  Soil  Response.  The  effect 
of  adjacent  track  segments  is  taken  into  account  by  the  factor  Ffl  which 
is  approximated  by  the  following  equation: 


Fa  =  Afa  +  Ffa  x  Dr 


where  Afa,  F^a  =  parameters 


The  factor  Fa  applies  to  the  calculation  of  the  dynamic  response  of  the 
ground  to  track  segment  "n"  when  either  of  the  adjacent  track  segments 
"n+1"  or  "n-1"  meets  the  following  conditions: 

•  The  adjacent  track  segment  is  in  the  direction  of  the 
interface  shear  stresses  acting  on  track  segment  "n" 
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•  The  position  of  the  adjacent  track  segment  is  lower  than 

that  of  track  segment  "n" 

•  The  velocity  of  the  adjacent  track  segment  is  downward. 

5.4.4.  Ground  Model  for  Use  with  TACOM  Dynamic  Track  Model.  The  ground 
model  consists  of  a  series  of  DYGREP  plates,’  as  shown  in  Figure  5-9. 
Realistic  simulation  requires  that  models  of  ground  response  be  fixed  to 
the  ground  in  their  horizontal  position,  as  opposed  to  ground  models 
which  travel  with  the  vehicle.  In  the  ground  model  consisting  of  a 
series  of  DYGREP  plates,  the  problems  arising  from  the  motion  of  the 
vehicle  relative  to  the  ground  model  are  resolved  by  adding  a  new  DYGREP 
plate  in  the  front  of  the  vehicle  and  omitting  one  at  the  rear  whenever 
the  distance  traveled  by  the  vehicle  equals  or  exceeds  the  width  of  a 
DYGREP  plate.  In  computer  methodology  this  is  achieved  by  an  overlay  of 
the  computer  model  of  DYGREP  plates.  The  vertical  position  of  omitted 
DYGREP  plates  may  be  preserved  in  the  program  if  it  is  desired  to  have  a 
record  of  the  rut  depth  made  by  the  tracks  of  the  vehicle. 

In  the  upper  part  of  Figure  5-9  a  potential  track  geometry  is  shown 
which  constitutes  (together  with  the  track  segment  velocities)  the. input 
to  the  ground  model  at  a  given  time.  This  input  geometry,  when  overlaid 
on  the  ground  model  of  DYGREP  plates,  shows  the  elevation  differences 
between  the  momentary  position  of  the  plates  and  track  segments. 

Whenever  the  elevation  of  a  track  segment  is  lower  than  that  of  the 
plate  beneath  It,  the  DYGREP  plate  is  forced  to  move  downward.  For 
computational  efficiency,  and  In  conformity  with  the  horizontal  rigidity 
of  the  DYGREP  plates  (which  does  not  allow  rotation  about  its  center), 
the  pin  displacements,  instead  of  those  of  the  track  segments,  have  been 
chosen  as  agents  which  activate  the  DYGREP  plates.  The  effect  of  pin 
displacement  is  assumed  to  extend  half  of  the  width  of  the  track  segment 
each  way  and  to  exert  a  displacement  constraint  on  the  ground  model  over 
this  range.  Generally,  a  pin  displacement  constraint  affects  two 
adjacent  DYGREP  plates,  and,  conversely,  the  forced  displacement  of  a 
DYGREP  plate  consists  of  a  proportional  allocation  of  two  pin 
displacement  constraints. 

The  sequence  of  calculations  for  the  determination  of  the  normal 
pressures  resisting  the  forced  displacement  of  a  DYGREP  plate  in  the 
ground  model  is  shown  in  Figure  5-10. 

The  tangential  (shear)  stresses  generated  at  the  track-soil  interface 
are  assumed  to  be  related  to  the  magnitude  of  slip  between  track  segment 
and  soil  by  the  following  equation: 

1  -  Traax  0  -  e'S/K)  l12' 

where  s  =  slip 

t  =  shear  stress  mobilized  at  the  track-soil 

interface 
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Dynamic  Soil  Response 


COORDINATES  OF  PIN  Mn" 

GROUND  COORDINATES  AT 
DVGREP  PLATE  T  (INITIAL) 


=  maximum  value  of  interface  shear  stress 


Tmax 

K  =  slip  parameter 

Slip  is  defined  as  follows: 


where 


X 

max 


*max  '  “  <*cg'  *tr> 


[13] 


Because  of  the  assumed  inextensibility  of  the  track  in  the  dynamic  track 
model,  the  velocity  of  the  individual  track  segments,  X,  ,  differs  but 
little  from  the  track  velocity.  Therefore,  for  simplicity  and  to 
minimize  computer  time,  the  track  velocity  is  used  for  the  determination 
of  slip. 


5.5.  Preprocessor  Program 

The  preprocessor  program  serves  the  purpose  of  computing  input  values 
for  the  dynamic  track-soil  interface  subroutine  from  experimental  data 
or  soil  properties  not  directly  usable  in  that  program.  The 
preprocessor  program  addresses  the  following  three  tasks: 

•  Determination  of  the  factors  Fs  and  Ffl 

•  Estimation  of  the  triaxial  strength  qt 

•  Determination  of  the  parameters  B  and  R  from  test  data. 

For  the  determination  of  the  factors  F$  and  F.  and  the  triaxial  strength 
qt,  the  track  segment  geometry  and  -the  strength  properties  of  the  soil 
(cohesion,  friction  angle  and  unit  weight)  are  needed  as  input  values  in 
the  preprocessor  program. 

5.5.1.  Determination  of  the  Factors  F.  and  Fa.  A  rectangle  of  the  size 
of  a  track  segment  is  subdivided  into  10  subrectangles  each  way  and  the 
bearing  stresses  at  the  nodal  points  (corners  of  subrectangles)  are 
computed  by  numerical  integration  of  the  differential  equations  of  the 
plasticity  theory  for  soils.  The  three-dimensional  case  of  a 
rectangular  plate  is  approximated  by  computing  bearing  stresses  for  two- 
dimensional  failure  in  the  x,-x,y,-y  directions  and  assuming  that  the 
three-dimensional  bearing  stress  is  the  lowest  one  of  the  four  computed 
for  a  nodal  point.  To  determine  the  relationship  between  the  Fs  factor 
and  the  interface  friction  angle  6,  the  three-dimensional  bearing 
stresses  are  determined  for  various  ratios  of  6  to  its  maximum  value 
6max*  was  f°und  that  a  parabolic  relationship  between  these  ratios 
ana  F$  approximates  the  computed  values  of  F$  reasonably  well.  The 
coefficients  of  the  best  fitting  parabolic  relationship  are  determined 
in  the  program  for  use  in  the  dynamic  track-soil  interface  model. 
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The  factor  Fa  is  determined  by  computing  the  three-dimensional  bearing 
stresses  for  the  same  ratios  of  6/6  as  for  the  computation  of  F$ 
but  for  a  rectangle  twice  the  width  8fxthe  track  segment.  It  was  found 
that  the  computed  Fa  values  for  the  various  6/6  ratios  can  be 
approximated  reasonably  well  by  a  linear  relationship.  The  coefficients 
of  the  best  fitting  linear  relationship  are  determined  in  the  program 
for  use  in  the  dynamic  track-soil  interface  model. 

5.5.2.  Estimation  of  Triaxial  Strength  q^.  The  triaxial  strength  of 
the  soil,  if  not  available,  may  be  estimated  from  the  soil  strength 
parameters  c  and  $  ,  and  the  unit  weight  of  soil  y  .The  chamber 
(lateral)  pressure  is  assumed  to  be  equal  to  the  lateral  pressure  in  the 
soil  at  a  depth  equal  to  half  the  width  of  a  track  segment.  The  soil  is 
assumed  to  be  loaded  at  the  surface  by  a  uniform  pressure  exerted  by  the 
weight  of  a  track  segment.  The  triaxial  pressure,  qt,  is  then  computed 
from  the  following  formula: 

qt  =  (o3  +  c. cotan  ($))  •  -  c.cotan($)  [14] 

where  c  =  cohesion 

41  =  friction  angle 
03  =  lateral  pressure 

5.5.3.  Determination  of  the  Parameters  B  and  R  from  Test  Data.  A  set 
of  test  data  pairs,  consisting  of  sinkage  values  in  inches  and  pressure 
values  in  psi,  is  the  input  for  this  program.  The  values  of  B  and  R  for 
this  data  set  are  determined  in  the  program  by  calculating  the  new 
variables  defined  in  Equation  [5]  and  determining  the  best  fitting 
values  by  linear  regression. 
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